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REDUCIBILITY OF MATRIX WEIGHTS 


JUAN TIRAO AND IGNACIO ZURRIAN 


Abstract. In this paper we discuss the notion of reducibility for matrix 
weights and introduce a real vector space Cm which encodes all information 
about the reducibility of W. In particular a weight W reduces if and only if 
there is a non-scalar matrix T such that TW = WT*. Also, we prove that 
reducibility can be studied by looking at the commutant of the monic orthog¬ 
onal polynomials or by looking at the coefficients of the corresponding three 
term recursion relation. A matrix weight may not be expressible as direct sum 
of irreducible weights, but it is always equivalent to a direct sum of irreducible 
weights. We also establish that the decompositions of two equivalent weights 
as sums of irreducible weights have the same number of terms and that, up to 
a permutation, they are equivalent. 

We consider the algebra of right-hand-side matrix differential operators 
T>(W) of a reducible weight W , giving its general structure. Finally, we make 
a change of emphasis by considering reducibility of polynomials, instead of 
reducibility of matrix weights. 


1. Introduction 

The subject of matrix valued orthogonal polynomials was introduced by M.G. 
Krein more than sixty years ago, see [211 [22]• In the scalar valued context, the 
orthogonal polynomial satisfying second order differential equations played a very 
important role in many areas of mathematics and its applications. In HU, A.J. 
Duran started the study of matrix weights whose Hermitian bilinear form admits 
a symmetric second order matrix differential operator, following similar considera¬ 
tions by S. Bochner in the scalar case in [2] . 

Motivated by previous works on the theory of matrix valued spherical functions, 
see m , one finds in [H HU HI] the first examples of classical pairs considered in 
HU. In 0 one can find the development of a method that leads to introduce some 
other examples of classical pairs. In the last years, one can see a growing number 
of papers devoted to different aspects of this question. For some of these recent 
papers, see (28] [9] 0] HU [24] [25] as well as HU HU E3 EUEU13 EU H] j among others. 

When one is in front of a concrete matrix weight one should ask himself: is 
it a new example or is it equivalent to a direct sum of known weights of smaller 
dimensions? In this work we prove that a weight W is reducible to weights of 
smaller size if and only if there exists a non-scalar constant matrix T such that 
TW = WT *, see Theorem 12.81 We also study the form in which a matrix weight 
reduces. In particular, we show that this reduction is unique. 

Given a matrix weight W and a sequence of matrix orthogonal polynomials 
{On}n£N 0 i it is natural to consider the matrix differential operators D such that 
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Q n is an eigenfunction of D for every n £ No- The set of these operators is a 
noncommutative *-algebra V(W) (see 5] OS [Ml 32]). In Section [3] we study the 
structure of the algebra T>(W) for a reducible weight, see Proposition 13.31 

Finally, we give some irreducibility criteria, showing that one can decide the 
irreducibility of a weight by also considering the commutant of the monic orthogonal 
polynomials {P„} ng N 0 (see Corollary 14.41) . or by looking at the commutant of the 
coefficients of the three term recursion relations satisfied by them (see Theorem 

USD. 

We prove in Theorem 14.121 that if there is a sequence of orthogonal polynomials 
with respect to an irreducible weight W and with respect to a weight V, both with 
bounded support, then we have IT = XV , with A > 0. This is not true when the 
support is not bounded, even more, one can have a sequence of polynomials which 
are orthogonal with respect to two weights W and V, being IT irreducible and V 
reducible (see Example 14.131) . Thus, at the end of the paper we make a change 
of emphasis by considering reducibility of polynomials, instead of reducibility of 
matrix weights. 

It is known that the classical scalar orthogonal polynomials do not satisfy linear 
differential equations of odd order. Even more, the algebras T>(W) are generated 
by one operator of order two, see [231 ■ The situation in the matrix context is 
very different, see [7 . Nevertheless, deciding whether there is or there is not an 
irreducible weight W, such that the algebra T>(W) contains an operator of first 
order, is a problem that is still unsolved. We believe that the results in the present 
paper may be useful for giving an answer in the future. 

Acknowledgement. This paper was partially supported by CONICET PIP 112- 
200801-01533 and by the programme Oberwolfach Leibniz Fellows 2015. 


2. Reducibility of Matrix Weights 


An inner product space will be a finite dimensional complex vector space V 
together with a specified inner product (•,•). If T is a linear operator on V the 
adjoint of T will be denoted T *. 

By a weight W = W ( x ) of linear operators on V we mean an integrable function 
on an interval (a,b) of the real line, such that W{x) is a (self-adjoint) positive 
semidefinite operator on V for all x £ ( a,b ), which is positive definite almost 
everywhere and has finite moments of all orders: for all n £ No we have 

[ x n W(x)dx £ End(V). 


More generally we could assume that IT is a Borel measure on the real line of 
linear operators on V, such that: W{X) £ End(T) is positive semidefinite for any 
Borel set A, IT has finite moments of any order, and IT is nondegenerate, that is 
for P in the polynomial ring End(T)[x] 


(P,P)= [ P(x) dW(x)P(x)* = 0 

J R 


only when P = 0. In any case the size of IT is the dimension of V. 

If we select an orthonormal basis {ei,...,ejv} in V, we may represent each 
linear operator IT( x ) on V by a square matrix IT (x) of size N , obtaining a weight 
matrix of size N on the real line. By this we mean a complex N x A-matrix valued 
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integrable function on the interval (a, b) such that W(x) is a (self-adjoint) positive 
semidefinite complex matrix for all x £ (a,b), which is positive definite almost 
everywhere and with finite moments of all orders: for all n £ No we have 

M n = f x n W{x) dx £ Matjv(C). 

J a 

Of course, any result involving weights of linear operators implies an analogous 
result for matrix weights and viceversa. Along the paper, we shall work either with 
weights of linear operators or with matrix weights, depending on what we consider 
better for reader’s understanding. 

Definition 2.1. Two weights W = W[x) and IT' = W'{x) of linear operators on 
V and V', respectively, defined on the same interval are equivalent , W ~ IT', if 
there exists an isomorphism T of V onto V such that W'(x) = TW(x)T* for all 
x. Moreover, if T is unitary then they are unitarily equivalent. 

Two matrix weights W = W{x) and IT' = IT'(x) of the same size and defined 
on the same interval are equivalent , IT ~ IT', if there exists a nonsingular matrix 
M such that W’(x) = MW(x)M* for all x. Moreover, if the matrix M is unitary 
then they are unitarily equivalent. 

Let S be an arbitrary set of elements. By an .S-module on a field K we mean a 
pair (IT, T) formed by a finite dimensional vector space V over I\ and a mapping 
W which assigns to every element x £ S a linear operator W{x) on V. It follows 
that an ^-module is an additive group with two domains of operators, one being 
the field K and the other one being the set S. 

In particular, if V is an inner product space and W = W(x) is a weight of 
linear operators on Id, then (IT, V) is an ^-module, where S is the support of W. 
It is important to notice that the equivalence of weights of linear operators does 
not coincide with the notion of isomorphism among S'-modules. If W = W(x) is 
a weight of linear operators on Id, we say that IT is the orthogonal direct sum 
IT = ITi ® • ■ ■ ® Wd if V is the orthogonal direct sum V = V\ ® • • • ® Tj where each 
subspace Vj is IT(ai)-invariant for all x and Wi{x) is the restriction of W{x) to V)- 

Definition 2.2. We say that a weight IT = IT(a;) of linear operators on V reduces 
to smaller-size weights, or simply reduces , if IT is equivalent to a direct sum IT / = 
W{ ® • • • ® W' N of orthogonal smaller dimensional weights. 

We say that a matrix weight IT = IT( x ) reduces to smaller-size if its correspond¬ 
ing weight IT = W(x) of linear operators reduces. 

Theorem 2.3. A matrix weight IT reduces to scalar weights if and only if there 
exists a positive definite matrix P such that for all x, y we have 

(1) W(x)PW(y) = W(y)PW(x). 

Proof. Suppose that IT(a:) = MA(x)M*, A = A(x) a diagonal weight. Then 

W(x)(M M*)~ 1 W(y) = (M A(x)M*)(MM*)~ 1 (M A(y)M*) = MA(x)A{y)M* 

= MA(y)A(x)M* = {MA(y)M*)(MM*)~ 1 (MA{x)M*) 

= W(y)(MM*)~ 1 W(x). 

Hence P = (M A/*) -1 is a positive definite matrix which satisfies (jTJ. 
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Conversely, assume CD- First of all take xq such that W(x o) is nonsingular. Let 
A be a symmetric positive definite matrix such that A 2 = W[x o). Take W'(x) = 
A~ 1 W(x)A~ 1 and P' = APA. Then W' is equivalent to W, P' is positive definite 
and W'(x)P'W(y) = W'(y)P'W'(x) for all x,y, and W'(xo) = I■ Thus without 
loss of generality we may assume from the start that W{x o) = I. 

Now the hypothesis implies that W(x)P = PW(x) for all x. Let E be any 
eigenspace of P. Then E is TF(a;)-invariant. Hence W(x) and W(y) restricted to 
E commute by m and are self-adjoint. Therefore W (x) restricted to E are simul¬ 
taneously diagonalizable for all x, through a unitary operator (cf. T9], Corollary 
to Theorem 30, p. 283). Since this happens for all eigenspaces of P and they are 
orthogonal we have proved that W(x) is unitarily equivalent to a diagonal weight, 
hence W reduces to scalar weights. This completes the proof of the theorem for a 
matrix weight. □ 

Theorem 2.4. The following conditions are equivalent: 

(i) A matrix weight W is unitarily equivalent to a diagonal weight; 

(ii) For all x,y we have W(x)W(y)=W(y)W(x); 

(in) There exists a positive definite matrix P such that for all x,y we have 
W(x)PW(y) = W{y)PW{x) and W{x)P = PW{x). 

Proof. 

(i) implies (ii) If U is unitary and W(x) = UA(x)U* where A(x) is diagonal for all 
x , then W{x)W(y) = (UA(x)U*)(UA(y)U*) = UA(x)A(y)U* = UA(y)A(x)U* = 
{U A{y)U*)(U A(x)U*) = W(y)W(x). 

(ii) implies (iii) It is obvious. 

(iii) implies (i) Let E be any eigenspace of P. Then E is TF(x)-invariant. Hence 

W(x) and W(y) restricted to E commute and are self-adjoint. Therefore W{x) re¬ 
stricted to E are simultaneously diagonalizable for all x, through a unitary operator 
(cf. [19], Corollary to Theorem 30, p. 283). Since this happens for all eigenspaces 
of P and they are orthogonal we have proved that W(x) is unitarily equivalent to 
a diagonal weight, hence W reduces to scalar weights. This completes the proof of 
the theorem. □ 

The following corollary can also be found in [51, p. 463]. 

Corollary 2.5. Let W = W(x) be a matrix weight such that W(xq) = I for some 
Xq in the support of W. Then W reduces to scalar weights if and only if W is 
unitarily equivalent to a diagonal weight. 


Proof. Let us consider first the case of a matrix weight. It is obvious from the 
definition that if W is unitarily equivalent to a diagonal weight, then W reduces 
to scalar weights. Conversely, let us assume that W reduces to scalar weights. By 
Theorem 12.31 there exists a positive definite matrix P such that W{x)PW(y) = 
W(y)PW(x) for all x, y. If we put y = x o we get W(x)P = PW(x) for all x. Then 
Theorem 12.41 implies that W is unitarily equivalent to a diagonal weight. □ 


Example 2.6. Let W{x) 



+ x 

X 


■A 

I supported in (0,1). Then 


W(x) 
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Therefore W reduces to scalar weights and W(x)PW(y) = W{y)PW{x) for all 
x, y G (0,1) with 


P = {MM*)- 1 



However W is not unitarily equivalent to a diagonal weight. In fact 

>™>=( x3 r i) 

({x 2 + x){y 2 +y) + xy (x 2 + 2x)y\ 

~ V (v 2 + 2 v) x 2x v )' 

Since the entries (1,2) and (2,1) are not symmetric in x and y, W{x)W{y) ^ 
W{y)W{x) for all x,y. 

Definition 2.7. We say that a matrix weight W is the direct sum of the matrix 
weights Wi, W 2 ,..., Wj, and write W = W\ © W 2 © • • • ® Wj, if 


W{x) 


(W\{x) 0 

0 W 2 (x) 


\ 0 0 


0 ^ 
0 


W 3 {x)j 


Theorem 2.8. Let W = W{x) be a matrix weight function with support S. Then 
the following conditions are equivalent: 

(i) W is equivalent to a direct sum of matrix weights of smaller size; 

(ii) there is an idempotent matrix Q ^ 0 ,1 such that QW{x) = W(x)Q* for all x; 
(in) the space Cr = {T G Matjv(C) : TW{x) = W{x)T* for a.e. x G S'} © K/. 


Proof. 

(i) implies (ii) Suppose that W{x) = MW(x)M* = W\{x) © W 2 {x) where W\ 
is of size k and W 2 (x) is of size n — k. Let P be the orthogonal projection onto 
the subspace of C” generated by {ei,...,efc}. Then PW(x) = W{x)P. Hence 
P(MW(x)M*) = {MW(x)M*)P and 


(. M~ l PM)W{x) = W (x)(M* P(M*)~ 1 ) = W(x)(M~ 1 PM)*. 


Therefore if we take Q = M 1 PM , then Q is idempotent, Q ^ 0,1 and for all x 
QW(x) = W(x)Q*. 

(ii) implies (iii) Clearly Cr is a real vector space such that / G Cr. Now the 
implication is obvious. 

(iii) implies (i) Our first observation is that if T G Cr then its eigenvalues are 
real. In fact, by changing W by an equivalent weight function we may assume that 
T is in Jordan canonical form. Thus T = © 1<i<s </?: is the direct sum of elementary 
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Jordan matrices Ji of size di with characteristic value Ai of the form 


/ A* 0 

1 Ai 


Ji = 


0 0\ 
0 0 


0 0 • • ■ Ai 0 

\0 0 • • • 1 Ai) 


Let us write W{x) as an s x s-matrix of blocks Wij(x) of di-rows and J,-columns. 
Then, by hypothesis, we have JiWu(x ) = Wu(x)J*. Thus, if 1 < k < n is the index 
corresponding to the first row of Wu(x), we have A iWkk{x) = Wkk{x)Ai. But there 
exists x such that W{x) is positive definite and so Wkk{x) > 0. Therefore Ai = Ai. 

Another important property of the real vector space Cr is the following one: if 
T £ Cr and p £ R[x] then p(T ) £ Cr. 

Let L and TV be, respectively, the semi-simple and the nilpotent parts of T. The 
minimal polynomial of T is (x — Ai) ri • • • (x — Aj) rj where Ai,. .., Aj are the different 
characteristic values of T and ri is the greatest dimension of the elementary Jordan 
matrices with characteristic value Ai. A careful look at the proof of Theorem 
8, Chapter 7 of [H] reveals that L and TV are real polynomials in T. Therefore 
L,N £ Cr. 

The minimal polynomial of L is p = {x — Ai) • ■ ■ (x — Aj). Let us consider the 
Lagrange polynomials 


Pk 


n 


{x - Aj) 

(Ai A k ) 


Since Pk(Ai) = dik, and L = L* it follows that Pi = Pi(L ) is the orthogonal 
projection of C” onto the Ai-eigenspace of L. Therefore Pi £ Cr. 

Let us now consider the nilpotent part TV of T: TV is the direct sum of the 
nilpotent parts TVi of each elementary Jordan block Ji, i.e. 


/0 0 
1 0 


Ni = 


0 0 
\0 0 


0 0\ 
0 0 


0 0 

f 0/ 


Since TV £ Cr we have TViWii(x) = Wu(x)N*. If TVi were a matrix of size larger 
than one, and if 1 < k < n is the index corresponding to the first row of Wu(x), we 
would have Wkk{x ) = 0 for all x, which is a contradiction. Therefore all elementary 
Jordan matrices of T are one dimensional, hence TV = 0 and T = L. 

If Cr D RI, then there exists T £ Cr with j > 1 different eigenvalues. Let 
M £ GL(TV, C) such that L = MTM~ l be a diagonal matrix and let Pi, 1 < i < j 
be the orthogonal projections onto the eigenspaces of L. Let W(x) = MW(x)M*. 
Then / = Pi + • • • + Pj, P r P s = PsPr = S rs P r ,l < r,s < j and P* = P r , 
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P r W{x ) = W{x)P r for all 1 < r < j. Therefore, 

W(x) = {P 1 + ---+P j )W(x)(P 1 + --- + P j )= PrW{x)P s 

= P l W{x)P 1 + • • • + PjW(x)Pj = mix) ® • • • © Wj(x), 

completing the proof that (iii) implies (i). Hence the theorem is proved. □ 

For future reference we introduce the following definition. 

Definition 2.9. If W = W(x) is a matrix weight function with support S , then 
we define the real vector space 

Cr = {T <E Mat at (C) : TW(x) = W{x)T* for a.e. ieS}. 

Any matrix weight W is equivalent to a matrix weight IF with the order zero 
moment M 0 = I. In fact, it suffices to take W = M ( f 1 ^ 2 IF Alff 1 ^ 2 . 

Proposition 2.10. Let W be a weight matrix with support S such that M 0 = I. 
Then its real vector space Cr is a real form of the commutant algebra 

C = {T £ Matjv(C) : TW{x) = W{x)T for a.e. x £ S'}. 

Proof. The hypothesis implies right away that T = T*, hence the proposition 
follows. □. 

It is worth to state the following immediate corollary of Theorem 12.81 which 
resembles a Schur’s lemma. 

Corollary 2.11. Let W = W{x) be a matrix weight . Then W is irreducible if and 
only if its real vector space Cr = R/. 

The following result can also be found in a recent work announced in [20]. 

Corollary 2.12. A matrix weight reduces if and only if the commutant of W = 
M 0 1 ^ 2 H r M 0 1 contains non-scalar matrices. 

Let ( W, V) be an abstract S'-module. Then (IF, V) is said to be simple if it 
is of positive dimension and if the only invariant subspaces of V are {0} and V. 
An S'-module is called semi-simple if it can be represented as a sum of simple 
submodules. 

Instead, when W = W(x) is a matrix weight (resp. a weight of linear operators 
on V) one says that W is irreducible when it is not equivalent to a direct sum 
of matrix (resp. an orthogonal direct sum of operator) weights of smaller size. 
Furthermore an operator weight will be called simple if has no proper invariant 
subspace. 

A semi-simple abstract S’-module (IF, V) can be represented as the direct sum 
V = Vi © • • • © Vj of a finite collection = {V)} of simple S-submodules. Moreover, 
if we have a representation of this kind, and if V' is any invariant subspace of V, 
then there exists a subcollection $ 0 of $ such that V is the direct sum of V' and 
of the sum of the submodules belonging to $ 0 - 

Conversely, let V be an S-module which has the following property: if V' is 
any invariant subspace of F, there exists an invariant subspace V" such that V = 
V' © V" . Then V is semi-simple (cf. Propositions 1,2 in Ch. VI, §1 of [6j). 

Proposition 2.13. Let IF = W(x) be a matrix weight with support S. Then the 
S-module (IF, V ) is semi-simple. 
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Proof. Let V' be any invariant subspace of V and let V" be its orthogonal 
complement. Then, since W(x) is self-adjoint for all x, V" is invariant. Therefore 
the .S-module (Wj V) is semi-simple. □ 

Let V = Vi © • • • ® Vj = V[ © • • • ® Vj, be two representations of an abstract semi¬ 
simple iS-module V as a direct sum of simple submodules. Then we have j = f 
and there exists a permutation er of the set {1,..., j } such that Vj: is isomorphic 
to V^j for all 1 < i < j (cf. Propositions 3, in Ch. VI, §1 of [6]). Clearly, this 
uniqueness result can be generalized to the following one: if V = V\ © • • • © Vj and 
V' = V{ ® • • • ® Vj, are two isomorphic S'-modules represented as direct sums of 
simple submodules, then we have j = j' and there exists a permutation er of the 
set {1,..., j} such that Vj is isomorphic to for all 1 < i < j. 

This last statement is not true for matrix weights or operator weights on an 
inner product space V, because, as we pointed out at the beginning, the equivalence 
among weights has a different meaning than the isomorphism of ^-modules. In fact, 
Example 12.61 illustrates this: the matrix weight 

"w-Cr;) 


has no invariant subspace of C™, i.e. it is simple, but it is equivalent to 


W'(x) 


xr 0 
0 x 


which is the direct sum of two scalar weights. Nevertheless the following holds. 


Theorem 2.14. Let W = IVj®- ■ • 8 W 7 and W' = W[®- ■ ■ 0VIC, be representations 
of the operator weights W and W' as orthogonal direct sums of simple (irreducible) 
weights. If W and W' are unitarily equivalent, then we have j = j' and there exists 
a permutation a of the set {1,..., j} such that IVj is unitarily equivalent to 
for all 1 < i < j. 


Proof. It is clear that it is enough to consider the case W = W'. We shall 
construct the permutation er. Suppose that a(i) is already defined for i < k (where 
1 < k < j) and has the following properties: a) a(i) ^ a(i') for i < i' < k; b) V( is 
unitarily equivalent to VL^ for i < k; c) we have the orthogonal direct sum 

V = ® K(i) © ® v i- 

i<k i>k 

We consider the invariant subspace 

£ = ®K W © ®^- 

i<k i>k 

Let E 1 - be the orthogonal complement of E which is, by Propositions 1 in Ch. VI, 
§1 of [6j mentioned above, the direct sum of a certain number of the spaces V(. On 
the other hand E ± is unitarily isomorphic to V/E, i.e to 14. It follows that E 1 - is 
simple and therefore E 1 - is one of the Vf, say E 1 - = V[ Q . Since C E for i < k, 
we have io ^ a(i ) for i < k. We define cr(fc) — io- It is clear that the function a(i) 
now defined for i < k + 1, satisfies conditions a), b), c) above, with k replaced by 
k + 1. 
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Because we can define the injective function cr on the set {1,..., j} we must have 
j' > j. Since the two decompositions play symmetric roles, we also have j > j', 
hence j = j'. The theorem is proved. □ 

If we come back to our Example 12.61 we realize that a matrix weight may not 
be expressible as direct sum of irreducible matrix weights. But such a weight is 
equivalent to one that is the direct sum of irreducible weights. Taking into account 
this fact and Theorem l2.8l we make the following definition. 

Definition 2.15. We say that a matrix weight (resp. an operator weight) is com¬ 
pletely reducible if it is equivalent to a direct sum of irreducible matrix weights 
(resp. an orthogonal direct sum of irreducible operators). 

Observe that Theorem 12.81 implies that every weight is completely reducible. 


Theorem 2.16. Let W = Wiffi- • -®Wj and W' = IT}®- ■ ■ 0 W-, be representations 
of the operator weights W and W' as orthogonal direct sums of irreducible weights. 
If W and W' are equivalent, then we have j = j' and there exists a permutation a 
of the set {1,..., j} such that Wi is equivalent to IT}^ for all 1 < i < j. 

Proof. Modulo unitary equivalence we may assume from the beginning that W 
and W' are matrix weights of the same size. Let W'{x) = MW ( x)M* for all x , with 
M a nonsingular matrix. We may write M = UP in a unique way with U unitary 
and P positive definite, and P = VDV* where V is unitary and D is a positive 
diagonal matrix. Then W'(x) = (UV)D(V*W(x)V)D(UV)*. Modulo unitary 
equivalences we may assume that W'{x) = DW(x)D. If we write D = D\®- ■ -®Dj 
where Di is a diagonal matrix block of the same size as the matrix block W t . 
Then W' = (DiW±Di) © • • • © ( DjWjDj). From the hypothesis we also have the 
representation of W' = W[ © • • • © JTj, as an orthogonal direct sum of irreducible 
weights. Now we are ready to apply Theorem l2.14l to conclude that j = j' and that 
there exists a permutation cr such that DiWiDi ss IT}^ for all 1 < i < j. Hence 
Wi ~ IW^.j for all 1 < i < j , completing the proof of the theorem. □ 

Theorem 2.17. Let W = W{x) be an operator weight equivalent to an orthogonal 
direct sum W\ © ■ ■ ■ © Wd of irreducible weights. If j(T) is the number of distinct 
eigenvalues of T £ Cr, then d = ma x{j(T) : T £ Cr}. Moreover, W is equivalent 
to a matrix weight W' = W[ © ■ ■ ■ © W' d where W[ is the restriction of W' to one 
of the eigenspaces of a diagonal matrix D £ Cr. Besides j(T) is the degree of the 
minimal polynomial of T. 

Proof. Suppose that W = W\ © • • • © Wd and let Pi,... ,Pd be the corresponding 
orthogonal projections. Define T = \\P\ + ■ ■ ■ + XdPd with Ai,..., Ad all different. 
Then clearly T S Cr. Hence d < max{j(T) : T £ Cr}. 

Conversely, let T £ Cr such that j{T) = max{j(T) : T £ Cr}. Modulo unitary 
equivalence we may assume that W is a matrix weight. In the proof of Theorem 
12.81 we established that T is semi-simple. Thus we may write T = A~ X DA with 
D a diagonal matrix. Let W'(x) = AW{x)A*. Then W' ~ W and D £ Cr(W / ). 
We may assume that D = D i © • • • © Dj^ T ^ where Di is the Jordan diagonal block 
corresponding to the eigenvalue Ai of D. Then W’ = W[ © • • • © IT'™ where the 
size of the block W[ is equal to the size of the block Di, for all 1 < i < j(T), 
because DW'(x) = W'(x)D. If some W[ were not irreducible we could replace 
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it, modulo equivalence, by a direct sum of matrix irreducible weights. Thus there 
exists a matrix weight W" ~ W' such that W" = W" ® ■ ■ • © W" is a direct sum 
of matrix irreducible weights with j(T) < j. 

By hypothesis W = W\ © • • • © Wd■ Hence Theorem 12.161 implies that d = j. 
Therefore d = j(T) and the W' are in fact irreducible. Moreover, there exists 
a permutation a of the set {l,...,d} such that W[ = W a ^y The theorem is 
proved. □ 


3. The algebra T>(W) of a reducible weight 

We come now to the notion of a differential operator with matrix coefficients 
acting on matrix valued polynomials, i.e. elements of Matjv(C)[x]. These operators 
could be made to act on our functions either on the left or on the right. One 
finds a discussion of these two actions in HU. The conclusion is that if one wants 
to have matrix weights that do not reduce to scalar weights and that have matrix 
polynomials as their eigenfunctions, one should settle for right-hand-side differential 
operators. We agree now to say that D given by 

D-'Em*), 3=4 

2—0 

acts on Q(x) by means of 

s 

QD = Y,d i (Q)(x)F i (x ) 

1=0 

Given a sequence of matrix orthogonal polynomials {Q n }neN 0 with respect to a 
weight matrix W = W(x), we introduced in [TS] the algebra D(W) of all right-hand 
side differential operators with matrix valued coefficients that have the polynomials 
Q w as their eigenfunctions. Thus 

(2) V{W) = {D : Q n D = T n (D)Q n , T n {D) G Matjv(C) for all n G N 0 }. 

The definition of V(W) depends only on the weight matrix W and not on the 
sequence {Q„}neN 0 - 

Remark 1. If W' ~ W, say W' = MWM *, then the map D i —> MDM~ X 
establishes an isomorphism between the algebras T>(W) and V{W'). In fact, if 
{Qn}neN 0 is a sequence of matrix orthogonal polynomials with respect to W, then 
{Q' n = M<5 ra M _1 } ng N 0 is a sequence of matrix orthogonal polynomials with respect 
to W'. Moreover, if Q n D = T n {D)Q n , then 

Q' n (MDM - 1 ) = {MT n {D)M~ l )Q' n . 

Hence MDM~ X G V(W') and T n (D') = MT n (D)M~ 1 . 

It is worth observing that each algebra V(W) is a subalgebra of the Weyl algebra 
D over Mat at (C) of all linear right-hand side ordinary differential operators with 
coefficients in Mat n (C) [x] : 

D = { D = &F '- : e Matjv(C)[®]}. 
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It is also interesting to introduce the subalgebra V of the Weyl algebra D defined 


by 


v = {d = Y, dip i e : de s^ < *}• 


A differential operator D is symmetric if (PD, Q) = (P, QD) for all P, Q G 
Matjv(C)[x]. The set <S(W) of all symmetric differential operators is a real form of 
V(W), i.e, 

(3) V(W) = S(W) © iS(W), 


see Corollary 4.5 of [IS]. 


In this section we will pay special attention to the algebra D(W) when the matrix 
weight W = W\ © W' 2 - Let N, Ni and N 2 be, respectively, the sizes of W, W\ and 
W 2 . If D G Djv we will write 


D = 


(Du 

\D 21 


d 12 \ 

D 22 ) ’ 


where D ^ is a right-hand size differential operator of size N, x Nj. 


Proposition 3.1. IfW\ and W 2 are equivalent weight matrices, then 
V(W X © W 2 ) ~ Mat 2 (X>(Wi)) ~ Mat 2 (C) © X>(Wi); 
both isomorphims are canonical. 


Proof. We first observe that by changing W by an equivalent weight matrix we can 
assume that W± = W 2 . Moreover, if {P n }„ eNo is the sequence of monic orthogonal 
polynomials with respect to W\ then {P n = P n ©P n }neN 0 is the sequence of monic 
orthogonal polynomials with respect to W. 

If D = ( f? 11 n 12 ) G 2?(W), then P n D = A n (D)P n , which in matrix notation 

\D 21 D 22 / 

is equivalent to 

(P n 0 \ (Du D 12 \ (AUD) A f 2 (D)\ (P n 0 \ 

Vo Pu)\d 2 1 d 22 ) V a 2! (D) A2 2 (d)J{o p n J- 

Therefore P n D^ = T” (D)P n for 1 < i,j < 2, which implies that D l? G D(W\). 
Conversely, given Dij G T>(Wi) for 1 < i, j < 2, then P n Dij = A n (D,j)P n . If we 

define T n = AJ^f) ’ then Pn ° = Hence ’ D e 

Thus, D G V(W) if and only if D y G V(W\) completing the proof of the 
theorem. □ 

Proposition [371] admits the following obvious generalization. 


Theorem 3.2. If W±, W 2 ,..., W n are equivalent weight matrices, then 
V(W 1 © • • • © W n ) ~ Mat„ (fDfWi)) ~ Mat„(C) © X>(Wi); 
both isomorphims are canonical. 


We move to consider the case of a direct sum of two arbitrary matrix weights W\ 
and W 2 not necessarily equivalent. Let {P ra } n eNo and {Q n }neN 0 be, respectively, 
the sequences of monic orthogonal polynomials corresponding to Wi and W 2 and let 
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Pn = Pn® Qn be the sequence of monic orthogonal polynomials of W = W\ © W 2 . 
If D £ V(W), then P n D = A n (D)P n , which is equivalent to 

(Pn 0 \ (Dll D 12 \ (A ft(P) A ? 2 (P)\ (Pn o \ 

^0 Qn)\D 2 1 D 22 ) V A 2 i(D) A% 2 (D))\0 Qn)- 

Therefore, P n Du = A r £ 1 (D)P n and Q n D 22 = A 22 {D)Q n which implies that Pn £ 
V{Wi) and that D 22 £ V(W 2 ). Besides, P„Pi 2 = A™ 2 (D)Q n and Q„P 2 i = 
A 21 (D)P n . Let 

DiVixMa = |-D = : Fi £ Mativ lX Ar 2 (C)[a;] j, 

i 

Pn 1 xn 2 = jP = £ Dat iX iv 2 : degP, < i|. 

i 

It is also interesting to introduce 

V(Wi, W 2 ) = {D £ Pat iX at 2 : P n D = A„(P)Q„, A n (D) £ MatAr lX jv 2 (C), n £ No}, 
T>(W 2l W\) = {D £ PATsxAy : Q n D = A„(P)P n , A n (D) £ Mat^xAr^C), n £ No}. 

When IT] = W 2 we have P(Wi, W 2 ) = V(W\). We also observe that T>(Wi, W 2 ) 
is a left ^(Wij-module and a right P(W 2 )-module, and that P(W 2 ,Wi) is a 
left P(W' 2 )-module and a right P(Wi)-module. Besides V(Wi, W 2 )V(W 2 , W\) C 
ViWPj and V{W 2 ,W 1 )V{W U W 2 ) C V{W 2 ). 

Clearly P £ P(W) if and only if Pn £ V{WP), P 22 £ T>(W 2 ), P 12 £ P(Wi, W 2 ) 
and P 21 £ P{W 2 , W\). In other words we have proved the following proposition. 


Proposition 3.3. IfW\ and W 2 are two matrix weights, then 


T>(Wi 0 W 2 ) 


( P(Wi) V(Wi,W 2 )\ 
[v(W 2 ,W 1 ) V(W 2 ) )■ 


It is worth to remark at this point that if one is interested in studying, for 
example, classical pairs {W, P) consisting of a weight matrix W and a second order 
symmetric differential operator P £ T> or the algebra V(W), it is not enough to 
consider irreducible weights. 

We also point out that the above proposition generalize to the case of a direct 
sum of r > 2 matrix weights, giving a similar structure for T>(W\ © • • ■ © W r ). 


3.1. An Example. 

The pair (W, D) considered here is borrowed from 5.1 in |15j by taking a = /3. 
Let W be the following matrix weight 

W{x) = *“(1 - *)“ x ~ 2x 1 + = x a (l - x) a (x Pi + P 0 ), 

supported in the closed interval [0,1] with a > — 1. Let us first compute the 
commuting space C = {T £ Mat 2 (C) : TW(x) = W(x)T* for all x} of W. A 2 x 2- 
matrix T = (Uj)ij £ C if and only if TF\ = F\T* and TFq = FqT *, or equivalently 
if and only if TI = IT* and TJ= JT* where / is the identity matrix and J is the 
anti identity matrix. This immediately leads to T = si + tJ , s,t £ R. Therefore, 
Theorem 12.81 says that W reduces to a direct sum of scalar weights and indicates 
how we can find them. 

A nontrivial idempotent in C is obtained precisely for s = 1/2 and t = ±1/2 Let 
us take T=(J + J)/2. The eigenvalues of a nontrivial idempotent transformation 
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are 1 and 0 and the corresponding eigenspaces are V\ = {Tv : v £ C 2 } and Vo = 
{v — Tv : v £ C 2 }. Hence the column vectors (1, l) 4 £ V\ and (1,-1)* £ Vo- If we 

^ ^ , then MTM~ X = diag(l,0). From the proof of Theorem 


take M 1 = , 

12.81 we know that W'{x) = MW(x)M* is a diagonal weight. In fact 


W\x) = \x a (l-x)« ( 
= x a (l-x) 


1 1 
1 -1 


— 2x T1 


—‘lx + 1 


1 1 
1 -1 


1— x 0 
0 x 


Let us consider the scalar weights Wi(x) = x a (\ — x)“ +1 and w 2 {x) = x a+1 {\ — x) a 
supported in [0,1]. Then W' = w\ © u> 2 - 

From 5.1 in |15j . by transposing, we take the following symmetric differential 
operator with respect to W: 

D = d 2 x(l -x) + d{X - xU) + V, 

"0 0 
,0 —2(a + 1)(1 + u ) j 

'2(a + 1) — u 0 

0 2(a T 2) -(- Uj 

u being a real parameter. Therefore D' = MDM~ X is symmetric with respect to 
W'. Then 

D' = d 2 x{l - x) + d{X' - xU') + V', 

'-i 1 


V = 
U = 


X = 


* + l-f 
-f-1 


U I 5 
2 / 


V' = MVM- 1 = (a + 1)(1 + u ) 
X' = MIAT 1 = 


1 -1 


a + 1 — 1 — u 

0 a + 2 


U' = MUM' 1 = 


2a + 3 — 1 — u 

— 1 — u 2a + 3 


If we write D' = I ^} l 


Dn 


21 


Dj 2 j, we have 


D' 1X = {d 2 x{l — x) + 9(a + 1 — (2 a + 3)ar) — (a + 1)(1 + u))Eu, 

D'i 2 = ( — 9(1 + u)(l — a;) + (a + 1)(1 + u))£t 2 ; 

Z? 2 i = (9(1 + u)x + (a + 1)(1 + u))F/ 2 ij 

£>22 = (9 2 x(l — x) + 9(a + 2 — (2a + 3)x) — (a + 1)(1 + u])E 2 2 , 

where {££,} is the canonical basis of Mat 2 (C). Let {p n }neN 0 and {<?n.}neNo be the 
sequences of monic orthogonal polynomials associated, respectively, to w i and w 2 - 
Then {P n = p n ® q n }n£ N 0 is the sequence of monic orthogonal polynomials with 
respect to W'. If we use Proposition 2.7 in [18] and P n D' = A ' n P n , then we get 

_ f—n(n + 2a + 2) — (a + 1)(1 + u) (1 + u)(n + a + 1) 

n (l + u)(n + a + l) —n(n + 2a + 2) — (a + l)(l + u)y 

In particular we obtain 

x(l — x)p" + (a + 1 — (2a + 3 )x)p' n + n(n + 2a + 2 )p n = 0. 
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As an immediate consequence of the hypergeometric representation of Jacobi poly¬ 
nomials 

n, n + a + /3 + 1 1 — y' 


P^Hv) = 2A 


a + 1 


we obtain 


Pn(x) = P n’ a+1 \ l ~ 

QnW = . [~o )n rL - 2*) 

(n + 2a + 2) n 

since the leading coefficient of is (n + a + /? + l)„/2 n n\. From the identity 

x) = (—l) n P^’ a (x) we get that 


q n (x) = (~l) n p n (l - x ). 


It is wise and rewarding to verify the identity: p n D 12 = A(D 12 )q n which is 
equivalent to 

-(1 - x)p' n {x) + (a + 1 )p n = {n + a + 1 )q n (x). 

By using the identity ^ Pn°= i(n + a + /? + l)P^“| 1 -/ 3 + 1 ^(a;) we see that (13.11) 
is equivalent to 

(n + 2a + 2)(1 - x)P^ 1 ’ a+2) (l - 2x) + (a + l)P.^ a+1 \l - 2x) 

= (n + a + 1)P^“ +1 ’ Q - ) (1 — 2x), 


or to 

n{n + 2 a + 2)(1 — x) 2 Pi 


—n +1, n + 2a + 3 
a + 2 


+ (a + l ) 2 2 Pi 


—n, n + 2 a + 2 
a + 1 


;x ) = (n + a + l) 2 2 Pi 


— n, n + 2 a + 2 
a + 2 


By equating the coefficients of P on both sides one can easily check that they are 
equal. 


It is worth to point out now the following facts: -D 12 P 21 G T>(w\) and 

AUPMi) = A' n {D' 12 )A' n {D' 21 ) = (1 + u y-(n + a + l) 2 . 

Also D' n € V{w\) and 

A^D'u) = — n(n + 2 a + 2 ) — (a + 1)(1 + u). 

If u 7 ^ —1, then A' n (D' n — (a + l)(a — u)) = —(1 + u)~ 2 A' n (D' 12 D 21 ). From here, 
and since the set of representations A^ [n € No) separates the elements of T>{w\) 
(Proposition 2.8, [18]), the following interesting factorization follows 

<9 2 x(l — x) + d(a + 1 — (2a + 3)x) — (a + l ) 2 = (9(1 — x) — a — l)(9x + a + 1). 

Similarly we have 

9 2 x(l — x) + 9(a + 2 — (2a + 3)x) — (a + l ) 2 = (dx + a + 1)(9(1 — x) — a — 1), 
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and these two second order differential operators are related by a Darboux trans¬ 
formation. When u ^ — 1 let 

-Du = ( d 2 x{l — x) + d(a + 1 — (2a + 3)x) — (a + l) 2 )i?ii £ V(w\), 

D \2 = ( — 9(1 — x) + a + l) E \2 £ V{wi,w 2 ), 

D 2 i = (dx + a + l) E 2 i £ T>(w 2 , w>i), 

D 2 2 = (p 2 x{\ — x) + d(a + 2 — (2a + 3)x) — (a + 1)”)£ , 22 £ D(^w 2 ). 


Proposition 3.4. T>(wi © W 2 ) can 6e presented as the complex algebra with gen¬ 
erators En, E 22 ,Di 2 , D 21 and the relations derived from the matrix multiplication 
table EijE rs = 6j r Ei S . Moreover 


T>{wi © w 2 ) = 


( V{wi) 

\p(w 2 ,Wi) 


V(wi,w 2 )\ 
V(w 2 ) J- 


If we define Du = D\ 2 D 2 i and D 22 = D 21 D 12 , then V[vji) is the polynomial 
algebra C[Du] andV(w 2 ) =C[I? 22 ]- While T>(wi,w 2 ) is the free left V(wi)-module 
(resp. the free right T>(pv 2 )-module generated by Di 2 . Similarly V(w 2 ,wi) is the 
free left 'D(w 2 )-module (resp. the free right D{wi)-module generated by D 2 i- The 
center of V(w lffiu^) is the polynomial algebra C[Du + D 22 ] ■ Moreover, T>(wi@w 2 ) 
is a free module over its center with basis {Eu, E 22 ,Di- 2 ,D 21 }■ 


Proof. In the classical cases of Bessel, Hermite, Jacobi and Laguerre scalar weights 
w = w(x) the algebra D(w) is well understood: it is the polynomial algebra in 
any second order differential operator £ T>(w) (cf. [24]). Therefore in our case 
V(wi) = C[L>n] and V{w 2 ) = €[£> 22 ]- 

Now we will show that T>(wi,w 2 ) has no differential operators of even order. In 
fact, if D £ T>(wi,w 2 ) were of even order, then DD 2 i £ D(w 2 ) would be of odd 
order. This is a contradiction which proves our assertion. Therefore we can assume 
that D £ T>(wi,w 2 ) is of the form 

2s+l i 

I> = 53 0VFJ (®), F i {x) = Y,x j F i j . 

2—0 j =0 

The coefficient of <9 2 ( s+1 ) in DD 2 i is xF 2s+ i(x) whose leading coefficient is F| s s i . 
But the coefficient of a differential operator in V{w 2 ) of order 2(s+l) is a polynomial 
of degree 2(s + 1). Therefore F^s+i 7 ^ 0- 

By induction on the order of a differential operator D £ D(wi,w 2 ) we will 
prove that V(wi,w 2 ) = T>(wi)Di 2 . We already know that in T>(wi,W 2 ) there is no 
nonzero constant differential operator. Let D £ V(wi,W 2 ) be of degree one. Then 
D = (d(xF( + F$) + Fq)E 12 with F(,F(,F° £ C. We consider D - F(D l2 £ 
T>(wi,w 2 )• Since the coefficient of d is constant it follows that D — -F 1 -D 12 = 0. 
Now assume that D £ V(wi,w 2 ) is of order 2s + 1 > 3. Then D = X^Xo 1 d l Fi(x) 
with F 2s +i{x) = (•t 2s+1 F 2 S S +i + * * *) Ei 2 , where the dots stand for lower degree 
terms and i^s+i G C x . 

We now consider D = D — (—l) s F^+i -D 12 £ T>(wi,w 2 ). By construction the 
degree of the polynomial coefficient of d 2s+1 in D is of degree less than 2s + 1. 
Therefore ord (D) < 2s + 1. By the inductive hypothesis there exists D' £ V(wi) 
such that D = D'Di 2 , completing the proof that T>(wi,w 2 ) = V{wi)Di 2 . Since it 
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is clear that D'D \ 2 = 0 implies £>' = 0 , we have proved that T>(wi,w 2 ) is a free 
left £>(wi)-module generated by D\ 2 . 

By observing that p(£>n)£>i 2 = p(£>i 2 £>2i)£>i2 = Di 2 p{D 22 ) for any p e C[x], 
it follows that D(wi,w 2 ) is the free right X>(u>2)-nrodule generated by D\ 2 . In a 
similar way can be establish that T>(w 2 ,wi) is the free left ©(-u^-module (resp. 
the free right X>(u>i)-module generated by D 2 \. 

Let D be an element in the center oiV{w\®w 2 ). If £> = pn(£>n)+pi2(£>n)£>i2+ 
P2i(D 22 )D 2 i +P22(D 22 ), Pu,Pi2,P2i,P22 G C[x], then 

DD 12 = Pii(Du)Di2 + P2i{D 2 2)D 2 2 G T>(wi,w 2 ) ® V(w 2 ), 

and 

D12D = D 12 P2i(D 22 )D2i + E12P22 (£ 22 ) G T>(wi) ® T>{w\, w 2 ). 

Hence p2i(D 22 ) = 0 and pn(£>ii)I?i2 = Di2P22(D 2 2 ) = p 2 2(Dii)Di 2 . Therefore 
p n = p 22 . Similarly we have, 

DD 2 i = Pi 2 (Dn)Dii + P22(D 2 2)D 2 i G T>(wi) © T>(w 2 , w 1 ), 

and 

D21D = D 2 ipii(Dn) + D 2 ipi2(Dn)Di2 G T>(w 2 ,wi) © V(w 2 )- 
We conclude that pi 2 (Dn) = 0 . Therefore, if the zero order term of D is zero, then 
D = pn(-Dn) +pn(£>22) =Pu(D n + D 22 ), because D n D 22 = D 22 Dn = 0 . Thus 
in any case D = p{Dn + D 22 ), P G C[x]. Conversely, if D = p{D n + D 22 ), then it 
is obvious that D commutes with all elements of T>(w\) and of T>(w 2 )- Moreover, 
DD12 = p(Dn + D 22 )Di2 = p(Dn)Di2 = Di2p(D 22 ) = Di 2 p{Dn + D 22 )- In the 
same way one proves that -D-D21 = D 2 iD. So, the center of V(wi © w 2 ) is equal to 
€[£>11 + £>22]- 

Finally we point out that, 

p(£>n + D 22 )Eu = p(£>n), p{Dn + D 2 2 )E 2 2 = p{D 22 ), 

P{Du + D 22 )Di2 = p(£>ll)£> 12 , p(D 11 + £> 22)£>21 = p(£> 22 )£> 21 - 

From here it is easy to see that {£n, £22, £>12, £>21} is a basis of V(wi © w 2 ) over 
the polynomial ring C[£>n + £>22]- The proposition is proved. □ 


4. Criteria 

In this final section we give some criteria to decide whether a weight is irreducible 
or not. First we state the following elementary lemma. 

Lemma 4.1. Let w = w(x) be an integrable function in a finite interval (a, b) such 
that x k w(x) are integrable for all k G No. If 

rb 


/ x k w(x) dx = 0, for all k G No, 
J a 


then w = 0. 


Proof. By hypothesis, if p G C[x] then f p(x)w(x) dx = 0. Given / £ C{[a, 6]) and 
e > 0, by Stone-Weirstrass theorem, there existsp G C[x] such that \f(x)—p(x)\ < e 
for all x G [a, b]. Thus, 

pb pb pb 

/ f(x)w(x)dx < / \f(x) — p(x)\\w(x)\dx < e / \w(x)\dx. 
fa J a J a 
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Therefore f^ f{x)w{x) dx = 0. The Borel measure defined by w is regular, hence 
by Riesz representation theorem w = 0. □ 


In the next theorem we give an alternative description of the real vector space 
Cr, introduced in 12.91 of a matrix weight W defined on a finite interval (a, b). 

Theorem 4.2. If W is a matrix weight of size N on a finite interval ( a,b ), let 
{-PnjneNo be the sequence of monic orthogonal polynomials and let Mq be the mo¬ 
ment of order zero ofW. Then 

C R = {T : TM 0 = M 0 T *, TP n = P n T for all n € N 0 }. 


Proof. If T £ Cr, then TW(x) = W(x)T* for a.e. x £ (a, b). Since Pq = I it is 
obvious that TPq = PqT. Now assume that n > 1 and that TPj[x) = Pj(x)T for 
all x and all 0 < j < n — 1. Then 

r b r b 


(4) {PnT, Pj) = J P n (x)TW{x)Pj (x)dx = J P n {x)W{x)T*P*{x)dx 


= {P n ,P 3 T) = (P n , TPj) = (P n , Pj)T* = 0. 


Hence P n {x)T = AP n {x ) for some A £ Mat at (C) and all x. Since P n is monic it 
follows that T = A. Therefore by induction on n we have proved that P n T = TP n 
for all n £ No- Moreover, it is clear that TMq = MqT*. Therefore 


Cr C {T : TM 0 = M 0 T*, TP n = P n T for all n £ N 0 }. 


To prove the reverse inclusion we start with a matrix T such that TM 0 = M 0 T* 
and TP n = P n T for all n £ No- We will first prove by induction on n £ No 
that TM n = M n T*. Thus assume that n > 1 and that TMj = MfT* for all 
0 < j < n — 1. If we write P n = x n I + x n ~ 1 A n _i + • • • + A 0 , then by hypothesis 
TAj = AjT for all 0 < j < n — 1. 

Now (P„, Po) = 0 is equivalent to M n + + • • • + AqMq = 0. Hence 

TM n = M n T*. Therefore 

f x n {TW{x) — W{x)T*)dx = 0 for all n £ No, 

J a 

which is equivalent to TW{x) = W{x)T* for almost every x £ (o,6), see Lemma 
ro Thus 

{T : TM 0 = M 0 T*, TP n = P n T for all n £ N 0 } C C R , 
completing the proof of the theorem. □ 


Theorem 4.3. If a matrix weight W supported on ( a,b ) reduces, then T>{W) con¬ 
tains non-scalar operators of order zero. If a, b £ R, then the converse is also 
true. 

Proof. Assume that W reduces, then there exists an inversible matrix M such 
that 

W{x)=M ^ W 2 {x)) M *' for every x e ( a ’ b )’ 
with W\, H '2 weights of smaller size. Then, it is clear that 

M (o o) M " 
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is an operator in T>(W). For the second statement, assume that a, b £ R and 
that there is a non-scalar order zero operator D £ T>(W). By ([3|), we have that 
D = T + iT', with T, T' £ V(W) symmetric differential operators. Since D is non¬ 
scalar then T or T' is non-scalar. Assume that T is non-scalar and let {P n }rieNo 
be the sequence of monic orthogonal polynomials with respect to W. Clearly, 
TP n = P n T for n £ N. Since the adjoint of T is given by T = ||P 0 || 2 T* ||P 0 || -2 
m Theorem 4.3]) we have that T £ Cr, see Theorem 14.21 If T is scalar, then T' 
is non-scalar. Analogously, it can be seen that T' £ Cr. Hence, Cr is non-trivial. 
Therefore, by Theorem 12.81 W reduces. □ 

Corollary 4.4. If a matrix weight W on (a, b) reduces, then the commutant of the 
monic orthogonal polynomials {P„} ng pj 0 , given by 

{T £ Mat n (C) : TP n (x) = P n (x)T, for all n,x}, 

contains non-scalar matrices. If a,b £ R, then the converse is also true. 

Theorem 4.5. Let W be a matrix weight on (a, b) and 

A n P n -i(x) + B n P n (x) + P n+ i(x) = xP n (x) 

be the three term recursion relation satisfied by the corresponding sequence of monic 
orthogonal polynomials {P n }neN 0 - IfW reduces, then there is a non-scalar matrix 
commuting with {A n , B n } ne t$ 0 . Even more, if a,b £ R, then the converse is also 
true. 

Proof. If W reduces then, by Corollarv l4.41 there is non-scalar matrix T commuting 
with P„ for every n, then we have 

A n TP n -\{x) + B n TP n {x) + TP n+ \(x) = T A n P n -i(x) + TB n P n {x) + TP n+ i(ij, 

which implies that, for every n, A n T = TA n and B n T = TB n . 

Assume a, b £ R and let T be a matrix commuting with A n and B n for all 
n £ No- Then, we have 

A n P n -i{x)T + B n P n (x)T + P n+1 (x)T = xP n (x)T , 

A n TP n -\{x) + B n TP n (x) + TP n+1 (x) = xTP n (x), 

for every n. Since Pq(x)T = TPq(x), by induction on n > 0 it is clear that 
P n (x)T = TP n {x) for all n £ No- By Corollary 14.41 the proof is finished. □ 

Corollary 4.6. Let {Qn}neN 0 a sequence of orthogonal polynomials with respect 
to a weight W and with respect to a weight V , both on a finite interval (a, b). Then, 
W is irreducible if and only if V is irreducible. 

We feel that it is time to introduce the following definition. 

Definition 4.7. Given a sequence of matrix polynomials {QnjneNo we sa Y that it 
is reducible if it is equivalent to a sequence {Qn}neN 0 (he. Q n {x) = MQ^xfM -1 , 
for all x,n and some non-singular constant matrix M) such that 

Qn{x) = ~ 0 ^ , for all n £ N 0 . 

V 0 Qn(x)J 

For more general definitions of similarity and equivalence of matrix polynomials 
the reader can see [12] . 
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Lemma 4.8. If {Q n }nG n 0 is a sequence of matrix orthogonal polynomials equivalent 
to {Q n = Qn © Q 2 }neN 0 > then {Q^jneNo an d {QnlneNo are sequences of matrix 
orthogonal polynomials. 


Proof. Let W be the orthogonality weight function of {Qn}neN 0 - Let M be a 
matrix such that Q n = MQ n M ~ x , for every n £ No. Thus, the sequence {Q ra }neN 0 
is orthogonal with respect to W = MWM*. Let 


W = 


(Wx,x 

V^2,l 


Wx, 2 \ 

W 2 J ’ 


it is clear that Wx,x and W 2 , 2 are also matrix weights. 
A direct calculation gives 


0 = (Qrn,Qn) 


! {QlrO Qn)w 12 

{QmiQn)w 2 i (QmiQn)w2,2. 


Therefore {Q^}neN 0 and {Q 2 }neN 0 are, respectively, orthogonal sequences with 
respect to the weight functions Wxp and W 2<2 - □ 


Theorem 4.9. A sequence of monic orthogonal polynomials {P n }raeN 0 with respect 
to a weight W is reducible if and only if its commutant, given by 

{T £ Mat n (C) : TP n {x) = P n (x)T , for all n,x}, 

contains non-scalar matrices. 


Proof. Necessity follows from Lemma 14.81 and Corollary 14.41 For the sufficiency 
we can assume, without loss of generality, that ||Po|| = I. We have that {P„}„ e N 0 
satisfy a three term recursion relation 


(5) A n P n _x{x) + B n P n (x) + P n+ x(x) = xP n (x), 

for n £ Nq. Let S n = ||P n || 2 . Then it is well known that flSJ) implies that 


( 6 ) 


S n =A n S n - 1 , B n S n = (B n S n )*. 


By hypothesis, there is a non-scalar matrix T commuting with P n for every 
n £ No- From the proof of Theorem 14.51 it is clear that T also commutes with 
A n and B n for all n £ No- Since So is the identity, it follows from (j6]) that T also 
commutes with S n for all n £ No- Therefore, every eigenspace of T is invariant by 
{Pm A n , B n , S n : n £ No}. Thus, we have that there exists a unitary matrix M 
such that, for every n £ No, 

= ( 0 I)' 

(A 1 4°\ /R 1 R°\ 

MA " M_1 = (o n 4J’ = ( o" 4J’ 

where P^{x),A J nl Bf, Sf are square matrices of size dj , for j = 1,2, and P°( x), A°, 
Bn,S® are dx x d 2 matrices. Since S n is positive definite we have that S X ,S^ are 
positive definite and 5° = 0. By ©, B n S n is Hermitian, thus P° = 0; also from 
(ED, it is clear that A° = 0. 

By conjugating ([ 5 ]) by M and looking at the upper-right blocks we see that 
KPu- l(*) + B l n P» + Pn+ x(x) = xP°{x), 
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for n £ No. Since Pq = I we have Pq = 0. Then, by induction on n we get 
P° = 0 for every n £ No- Thus MP„M _1 = P* ® P 2 , completing the proof of the 
theorem. □ 

Corollary 4.10. A sequence of orthogonal polynomials {Qn} n eN 0 with inspect to 
a weight W on (a, b ) is orthogonal with respect to a reducible weight if and only if 
there is a non-scalar operator of order zero in D(W). 

Lemma 4.11. Let {Qn}ragN 0 be a sequence of orthogonal polynomials with respect 
to a matrix weight W and with respect to a matrix weight V. If 

f W(x)dx = I V(x)dx, 


then W and V have the same moments. If both weights have compact support, then 
W = V. 

Proof. Given the sequence {(3n}neN 0 : let us consider the corresponding sequence 
of rnonic orthogonal polynomials {P n }neN 0 - This sequence is also orthogonal with 
respect to W and V, therefore they satisfy a three term recursion relation 

A n Pn— 1(*^) T d3 n P n {x) -\- P n+1 (x) — xP n {pc) 


for n £ Nq. From the hypothesis we have that ||Po||^ = ||Po||y- Besides, from 


1 2 w = \\Po\\l 

we have ||P„+i|| 2 = A„ + i||P„|| 2 for every n £ No and for both norms 
|| • ||Therefore ||P„||^ = ||P„||y for any n £ N 0 . Hence, 


I w and 


P n (x)(W{x) - V {x))Pf n {x)dx = 0, 


for every n, m > 0. Hence, W and V have the same moments. If both have compact 
support, then by Lemma ri.il we have W = V. □ 

Theorem 4.12. Let {Q„}ngN 0 be a sequence of orthogonal polynomials with respect 
to an irreducible weight W and with respect to a weight V, both with compact 
support. Then, V is a scalar multiple ofW. 

Proof. Let us consider the corresponding sequence of monic orthogonal polynomials 
{P n }n g N 0 . From © we have that ||P„+i|| 2 = A n+ \ ||P n || 2 , for n £ No and the inner 
products given by W or V. As a straightforward consequence we obtain 


■MPo\\ 2 y = 


\Pn\\ 2 V, 


An ■ ■ ■ Ai ||Po||w — l|-fn||vy> -^n 
for £ No- Hence, since the norms are nonsingular hermitian matrices, we have 
II Pnfw IIPoll^ 2 = \\Pn\\ 2 v \\Po\\v 2 and ||P 0 ||^ 2 ||P„|| 2 W = ||P 0 ||^ 2 \\P n fv 
Then, 

An... AlllPoll^ = HPofy llPollv 2 \\Pn\\ 2 v = \\Po\\ 2 V \\Po\\w \\Pn\\w 

= ||Po||^||Po||^A„...A 1 ||P 0 || 2 f . 

Hence, ||Po||y ll-Pollw commutes with A n ■ ■ ■ A\ for every n £ No- Since A n is 
nonsingular for all n > 1 it follows that ||Po||y ||Po||^T commutes with A n for all 
n £ Nq. 


REDUCTIBILITY OF MATRIX WEIGHTS 


21 


For n € No and the inner products given by W or V we have Pn||Pn|| 2 = 
(B n \\P n \\ 2 )*, see ©. Thus 

BJPnfw = {Bn\\P n \\wY = (- Bn\\Pn\\ 2 V \\P 0 \\y 2 ||P 0 ||^)* 

= \\Po\\ 2 w \\p 0 \\v 2 (B n \\p n r v y = \\p 0 r w \\p 0 \\-*B n \\p n r v 
= l| J Po||^||Po||^B„||Po||^||Po||^||Pn||^. 

Therefore, ||Polly|| Polly/ 2 also commutes with B n for every n G No- By assumption 
W is irreducible, thus Theorem 14.51 implies that ||Po||y = A||P 0 ||^, with A> 0. 
Finally, Lemma Pi .III completes the proof of the theorem. □ 

When the support is bounded, we have proved that the reducibility of a weight 
is equivalent to the reducibility of any of its sequence of orthogonal polynomials. 
Even more, we showed that there is only one, up to scalar, irreducible weight for 
a given irreducible sequence of orthogonal polynomials. When the support is not 
bounded, this is not true; this is one of the main motivations behind Definition Id. 71 
We give the following ad hoc example to help the reader to understand this. 


Example 4.13. Let us consider the weight w defined by 

w{x) = —e~( los ^ / 2 , for x > 0, 
x 

and let {pn}neN 0 & e a sequence of orthogonal polynomials with respect to w. Let 

/( x) = ie - ( losx ) / 2 sin( 27 r log a;), for x > 0 . 

The function x k f(x) is continuous in ( 0 ,oo) and \x k f(x)\ < ^-e~^ ogx ^ 2 . By 
making the change of variables log x = s + k we get, 


' x k e -(\ ogx) 2 /2^ = e k 2 /2 


- s2 / 2 ds = 




0 fe 2 /2 


Hence x k f(x) is an integrable function in (0, oo), for all k > 0. Moreover, 

/* oo n oo j 

/ x k f{x)dx= / x k e~( logx ' ) / 2 sin( 27 rloga;) — 

Jo Jo x 


= e fe2 / 2 


o-P/2 


sin( 27 rs) ds = 0 . 


Now, we define the matrix weight W by 

r 2w(x) f(x) 


W(x) = 


/(a 


j{x) 


for x > 0. 


The reader can easily check that W(x) is positive definite for all x > 0. Also, 
by Theorem l£.<Sl after some straight forward computations, it follows that W is 
irreducible. Then, the sequence of polynomials {Pn}neN 0 > given by 


P n {x) 


(Pn(x) 0 A 
V 0 Pn(x))’ 


is orthogonal with respect to the irreducible weight W. Clearly, {PnlngNo °J S0 
orthogonal with respect to the matrix weight 


V(x) 


(w(x) 0 \ 

y 0 w(x)) ’ 


for x > 0. 
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